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Abstract 



We show that the axion decay constant does not get any correction at any 
order by external magnetic fields. On the other hand, in the context of the 
Wilsonian effective action under external magnetic fields, the axial currents 
get a finite correction. We then calculate the effect of strong magnetic fields 
B > 10 18 C7 on the axion- nucleon coupling and find that if B > 10 20 G in 
strongly magnetized neutron and white dwarf stars, the emission rate of axions 
is enhanced by several orders of magnitude. 
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The stellar evolution tightly constrains the emission of weakly interacting, low-mass 
particles such as neutrinos or axions [|l| whose properties are otherwise hard to probe ||. 
Recently, the emission rate of such particles from stars has been calculated more accurately 
by including the many-body effects ||, the effect of strong magnetic fields [|||, and the 
thermal pions and kaons ||[7|. The many-body effects in the dense core of stars are found 
to suppress the axion or neutrino emission rate by a factor of about 2, while the thermal 
pions and kaons to increase the axion emission rate by a factor of 3-4 within a perturbative 
treatment. Under a strong magnetic field B the cyclotron emission of neutrinos and axions 
offers a new cooling mechanism for magnetized white dwarfs and neutron stars. It puts 
a bound g ae < 10 -10 on the axion electron coupling constant. Similarly from magnetic 
white dwarf (WD) stars g ae < 9 x 10~ 13 (T/10 7 K) 5/4 ( 5 /10 10 G)~ 2 , where T is the internal 
temperature of the white dwarf. 

In this letter, we calculate the corrections to the emission rate of axions by nucleons in a 
strongly magnetized neutron or WD star whose magnetic field is stronger than the critical 
field, B > B c , such that the Landau gap is larger than the rest mass energy or Aq C d- As 
we will see later, for B > 10 20 G (or \/\qB\ > 4Aqcd — 800 MeV) the axion emission is 
enhanced quite a lot and gives a stringent bound on the axion-nucleon coupling. Large 
magnetic fields B ~ 10 12 -10 14 G have been estimated at the surface of neutron stars from 
the synchrotron radiation spectrum || and fields B ~ 10 18 -10 20 G are predicted in the core 
by scalar virial theorem ||. 

Let us first consider how the axion decay constant is renormalized in QED under a strong 
magnetic field. The calculation will go in parallel in QCD also. Under a constant magnetic 
field oriented along the z axis B = Bz, the spectrum of charged particles of charge q and 
mass m is quantized by the Landau level 

E A = a^Jm 2 + k 2 z + 2\qB\n, (1) 

where a = ± is the sign of the energy and the quantum number n labeling the Landau level 

2n = 2n r + 1 + \m L \ - sgn(qB)(m L + 2(3) (2) 

with /3 = ±2 the spin component along the magnetic field. Here n r is the number of 
nodes of radial eigenfunction and mi is the angular momentum. (Note that there will be 
corrections to the energy eigenvalue Eq.([l|) due to QED loop effects like the anomalous 
magnetic moment, but they are irrelevant since the corrections are smaller than the Landau 
gap.) 

At energy lower than the Landau gap, E < ^J\qB\, the relevant degrees of freedom in 
QED are photons and the lowest Landau level (LLL) fermions. Their interaction can be 
described by an effective action derived by integrating out the modes at the higher Landau 
levels (n > 0) [ID]. The effective action contains a (marginal) four-fermion interaction, which 



leads to chiral symmetry breaking even at weak attraction [I I - [131 . Using this effective 
action, we calculate the axion decay amplitude to two photons. The Wilsonian effective 
action of QED is given as [K 
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where the coefficients a\,b\, and c\ are of order e 2 while g\ is of order e 4 , determined by the 
one-loop matching conditions, and the ellipses are the higher order operators in momentum 
expansion. The components d± (d\\) are perpendicular (parallel) to the magnetic field and 

~ M = f i7"7V In 2 sgn( eJ B) if // = 0, 3 , , 

7 if// = 1,2. 1 ' 

Now, we consider the axion coupling to electrons: 

Ant = 7T^— d"a*7„76*, (5) 

where a is the axion field and /pq is the axion decay constant. We integrate out all the 
modes (^ n ^o) except the LLL electron (■0) to get the low energy effective interaction for a 
and ip, which is given as 

4£ = ^-PaftysKl, - 2 \ e B\ fpQ ^ a ^d a a^i d, /Mj + • • • , (6) 

where the ellipses denote the effective interaction terms containing the higher powers of 
axions and momenta, generated by the exchange interaction of \l/ n ^o- 

In the leading order, there are three diagrams that contribute to the axion decay am- 
plitude (see Fig. 1). We find that the first two diagrams vanish but the third diagram 
gives 

ie 2 f r ~ ~ i 

Ai = -j^j-(2n) 4 6\k -P)J(q + k) u k a e,( Pl )e^p 2 )Ti [-fS(q + Pih^sf 7 /3 %)J (7) 
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8\k - P)e^ Plfl e v ip 1 ) P2a e f3 {p 2 ), (8) 



where P = p± + P2 and e /Lt (p) the photon wave function, and the LLL electron propagator 

S(l) = -^—P_ e - l V\ eB \ (9) 
/II ~ m 

with P_ = 1 — i7 1 7 2 sgn(e5). The effective Lagrangian for the axion decay amplitude is then 

2 

£ eff = — —FF (10) 

We find that the axion decay constant is not renormalized in the leading order under an 
external magnetic field, fpQ(B) = fpq(0), in contrast with the pion decay constant that 
increases under an external magnetic field [|14]]. In fact, as is shown in the appendix A, 
the non-renormalization of the axion decay constant under external magnetic fields is exact 
simply because the axion couples to E ■ B not to magnetic fields B alone. 

Similarly, we calculate the decay amplitude of the axion into one photon (Fig. 2), using 
the effective Lagrangians Eqs. (||) and @, to get 
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l JPQ Jq 
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(2n)'S A (k- P y 2 ^e a (p)k B. 



(11) 
(12) 



^ 2 fp Q 

This result agrees with the direct calculation done with the effective Lagrangian Eq. (|lOl), 
where one of two photons is replaced by the external photon, assuring that fpg is not 
renormalized. We see that the external magnetic field does not change the axion-photon 
interaction. Therefore, the energy-loss processes involving only axions and photons such as 
the axion-photon conversion or the Primakoff process do not get affected by the dimensional 
reduction. 

Now we consider the process that electrons or quarks emit axions. For an extremely 
strong magnetic field B > B c , all the charged particles are in the lowest Landau level and 
thus the Bremsstrahlung emission is more relevant than the cyclotron emission, unless the 
interior temperature of stars is comparable to the Landau gap. The one loop correction to 
the axion Bremsstrahlung amplitude is (see Fig. 3) 
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(13) 
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We see that the Bremsstrahlung process has an extra contribution coming from the effective 
four-fermion interaction. (Note that here we assume k 2 <C m 2 . But, if m — > 0, the contribu- 
tion goes to zero.) In QED, g± oc e 4 and therefore quite small [10|, but in QCD, for a strong 
magnetic field B > 10 19 G (or \J\qB\ > Aqcd) such that the effective action for QCD is 
a 2 and is no longer negligible (see the appendix B for details). Especially, for 



valid, 



9x 



> 



quarks under a strong external magnetic field B > 10 20 G, g{ will be quite large. Thus 
the four-quark interaction will be stronger than the one-gluon exchange interaction and will 
break the chiral symmetry of QCD at energy above Aqcd- 

The axion emission process of neutron stars was first calculated by Iwamoto and it is 
found that axion bremsstrahlung from nucleon-nucleon collision is the dominant energy- 
loss process in neutron stars |TS]. Since the emission rate is proportional to the square of 



the axion-nucleon coupling g 2 n , any correction to this coupling will affect the emission rate. 
What we calculated above can be thought of as a finite correction to the axial singlet current 
of quarks and equally valid to the axial octet quark currents. Namely, both g7 M T a 7 5 g and 
qi^l^q have the same correction (1 + g{/ (27r 2 )). As was discussed in Ref. the derivative 
couplings of the axion to the axial vector baryon current is 



j— |2tr (Q A T a ) (FtrS 7 " 75 [T a , B] + DtvB^ l5 {T a , B}) + -trQ A tr J B 7 ^7 5 5 



(15) 



where B is the baryon octet matrices and Qa the axion charge. According to the current 
algebra, F, D, and S are defined to be proportional to the nucleon matrix elements of the 



axial quark currents [17]. Therefore we see that the axion-nucleon coupling gets a correction 
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2tt 2 ' 



(16) 
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Besides the correction due to the four-fermion interaction, the emission rate gets correc- 
tions since the axion-nucleon coupling and the nucleon mass will change due to the strong 
magnetic field. As a good approximation, we can take the nucleon mass as 



m n oc (qq} 1/3 . 

In chiral perturbation theory, the quark condensate is given as 

b dS V3iC (m q , B) 



(qq) 



m a =0 



(17) 



(18) 



where £ vac is the vacuum energy density. The one-loop vacuum energy under a magnetic 
field was calculated by Schwinger |L9| as 



^vac(-B) = £vac(0) — 



167T 2 



ds 



-Mis 



eBs 



sinh(e-Bs) 



(19) 



Therefore, the condensate under an external magnetic field is, using the Gell-Mann-Oakes- 
Renner relation F%M% = (qq) (m u + m d ), 



(qq) B = ( qq ) B 



1 + 



|efi|ln2 
16tt 2 F 2 



+ 



(20) 



where the ellipses denote the higher loop corrections |TJ|. Since the axion-nucleon coupling 
is proportional to the nucleon mass (g an oc m n j /pq) and the axion decay constant does not 
change under magnetic fields, 



9an{B) = g an (0) 1 + 



|eS|ln2 



(21) 
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The energy-loss rate per unit volume due to the nucleon-nucleon axion bremsstrahlung is 



calculated in Ref. [15 



and can be written as 



J_ 



2.5 2 rrQ.5 
m n 9an T 



(22) 



where / is the pion-nucleon coupling and T is the temperature. Therefore, we find the 
correction due to the strong magnetic field is 



Ql(B) = Qi(0) 1 



9l_ 
2vr 2 



+ 



1 + 



\eB\ In 2 
16tt 2 F 2 



+ 



2.5/3 



(23) 



where we have used the Goldberg- Treiman relation m n = fF n together with the Gell-Mann- 
Oakes-Renner relation. Similarly, the lowest-order energy emission rate per unit volume by 
the pion-axion conversion tt~ + p — >• n + a |]] has corrections due to the strong magnetic 
field as 



q: (B) = q: ( o)(i + |^+. 



|efi|ln2 
16tt 2 F 2 



+ ■ 



-1/3 



(24) 
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Note that the above derivation of the change in the quark condensate relies on the chiral 
perturbation theory and the one loop vacuum energy calculation, which will not be valid 
presumably for an extremely strong magnetic field such that the Landau gap becomes bigger 
than the p meson mass v/jei?j > m p . But, the condensate will still increase as B even for 
such a strong magnetic field because the increment of vacuum energy due to the external 
field will be balanced by forming larger condensates. 

The correction to the emission rate due to the external magnetic field comes from two 
sources. First one is from the finite correction to the axial current due to the effective 
four-quark operator and the second one is by the change of the quark condensate under 
the external field. For B = 10 18 ~ 10 19 G, the correction by the second one is less than 
a few percent while the correction by the first one is of order one. For B = 10 20 G or 
V |eB| = 0.77GeV, the correction by the second is about 24% for the Bremsstrahlung 
and about —10% for the pion-axion conversion, while the correction due to the first one 
is about 0.06a s for both processes, because the four-quark operator will evolve sufficiently. 
For B > 10 20 G the perturbative calculation breaks down and we expect that the correction 
will be extremely large because the strong enhancement of g{ at low energy. Therefore, if 
B > 10 20 G in stars, they will emit axions too quickly and will not survive present. On the 
other hand, if we observe such a strongly magnetized star, the axion-nucleon coupling has 
to be smaller than the current upper bound by several orders of magnitude. 

In conclusion, we have shown that the axion decay constant does not get a correction 
by the external magnetic field and we have calculated the correction to the axion-nucleon 
coupling due to external magnetic fields. We find that if B > 10 20 G in strongly magnetized 
neutron and white dwarf stars, the axion emission rate is enhanced by several orders of 
magnitude. 
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APPENDIX A: QED VACUUM ENERGY 



The low energy effective action of QED, integrating out fermions and scalars, has been 
calculated for various external fields [p2]| . Here, we calculate the QED vacuum energy, 
integrating out the axion fields: 



e Jx 



J x S v , c (E,B) 



d[a] exp 



1.2 l, Q _x2 



-F - - (d^ay + V cos 
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f PQ J + 16vr 2 f PQ j 



ie 2 aFF\ 



(Al) 



where we neglect other matter fields for simplicity but it does not change our argument. 
The vacuum energy density is then 
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£vac{E, B) = ^{E 2 + B 2 ) + E*JE, B), (A2) 

where the second term is due to axions and obviously a function of FF or E ■ B. If E ■ B = 
0, the vacuum energy does not change and neither does the axion decay constant. But, 
if E ■ B 7^ 0, the vacuum energy increases by the axion loop contribution because the 
axion-photon coupling in Euclidean space is imaginary. Using Schwarz inequality, one can 
easily show that 8^ 3JC (E,B) > 0. Thus under the external fields the condensate has to 
increase to balance the vacuum energy and the axion decay constant gets bigger. For general 
electromagnetic fields, we get f PQ (E,B) = f PQ {E ■ B) > f PQ {0). 

If we keep only the axion mass term in the potential in Eq. ( TAT| ) and integrate over a, 
we get for slowly varying fields 

OT ^) 2 (|£) 2 , (A3) 

where we used the axion mass m a = Aq CD / Jpq. 



APPENDIX B: QCD UNDER A CONSTANT EXTERNAL MAGNETIC FIELD 



Since quarks have electric charges, magnetic catalysis will operate for quarks as well when 
the external mag netic field is sufficiently strong, B > 10 19 G or y/\qB~\ > A QCD ~ 200MeV. 
By the renormalization group analysis, we will see that if the magnetic field is stronger than 
10 20 G, the four-quark interaction is stronger than the one-gluon exchange interaction and 
thus the chiral symmetry will break at scale higher than Aqcd • The derivation of low energy 
effective Lagrangian will be same as that of QED except the group theoretical factors. At 
low (but quite higher than Aqcd) energy the relevant degrees of freedom are LLL quarks, 
gluons, and photons. But, since the gluons couple to quarks more strongly than photons, 
we will neglect the photon interactions. As in QED, the exchange of quarks at the higher 
Landau level will generate a tree-level interaction of quarks and gluons, which will induce a 
four-quark interaction at one- loop matching (see Fig. 4); 



Lcs3 4\qB\ 



(QoQo) 2 + (Q0275Q0) 2 , (Bl) 



where Qo is the LLL quark and q is the electric charge of quark. The one-loop RG equation 
for the four-quark interaction is given as 

A*^ = ~ oJ(1b2) 2 , (B2) 

where a s is the coupling constant of strong interaction. Unlike QED, the strong coupling 
constant runs at scale below the Landau gap A^ = y/\qB\ as pi] 

1 1 11 , 

—r -r = /A r + — In -p. B3 

a 8 (fj.) a s (A L ) 2vr A L 
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Note that quarks do not contribute to the running coupling since the quark loop is finite 
due to the dimensional reduction. Combining the above RG equations (|B2|) and (B3), we 
find how the four-quark coupling changes as scale; 



</£(//) = 1.1424 (a s (/i) - a s (A L )) + g{{A L ), (B4) 

where gl(A L ) is of the order of a s (A L ) 2 . We see that the four-quark coupling becomes 
stronger than the strong coupling constant at a scale larger than the QCD scale fi > Aqcd 
if Al > 4Aqcd or B > 10 20 G. For such a strong magnetic field, the chiral symmetry of 
QCD will be broken at scale larger than Aqcd, the usual chiral symmetry breaking scale of 
QCD. 
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FIGURES 





FIG. 2. Axion conversion into a photon. Solid line denote electrons. 
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FIG. 3. Axion Bremsstrahlung emission. Solid lines denote electrons or quarks 
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(b) 

FIG. 4. Matching conditions. Double lines denote quarks at the higher Landau levels (n > 0), 
single line LLL quarks, and curled lines gluons. (a) tree level, (b) one loop level 
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